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In the first part of this work the classical and statistical aspects of the 
dynamics of an inextensible chain in three dimensions are investigated. In 
the second part the special case of a chain admitting only fixed angles with 
respect to the z— axis is studied using a path integral approach. It is shown 
that it is possible to reduce this problem to a two-dimensional case, in 
a way which is similar to the reduction of the statistical mechanics of a 
directed polymer to the random walk of a two-dimensional particle. 

PACS numbers: 05.40.-a, ll.10.Lm, 61.25.H- 

1. Introduction 

In this report we investigate the dynamics of an inextensible three- 
dimensional chain fluctuating in some medium at fixed temperature T. The 
chain is considered as the continuous limit of a freely jointed chain, which 
consists of a set of N — 1 rigid links of length a and N beads of mass m 
attached at the joints between two consecutive segments. The formulation 
of the dynamics of a chain with rigid constraints based on the stochastic 
equation of Langevin has been extensively studied in a series of seminal pa- 
pers by Edwards and Goodyear [H |3] . Unfortunately, to deal with these 
constraints at the level of stochastic equations is a cumbersome task. Up 
to recent times, most of the developments in the dynamics of a chain with 
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rigid constraints have been confined to numerical simulations, see for exam- 
ple Refs. [H|5j[6]. For this reason, we recently proposed an interdisciplinary 
treatment to this problem combining methods of field theory and statistical 
mechanics [7, . The strategy is to regard the change of the chain conforma- 
tion as the motion of Brownian particles with constrained trajectories. The 
framework of the calculations is that of path integrals. The constraints are 
introduced by a procedure which is commonly applied in statistical mechan- 
ics in order to enforce topological conditions on a system of linked polymers. 
One ends up in this way with a field theory which is a generalized non-linear 
sigma model (GNLcM). Recently, this path integral formulation has been 
connected to the usual description of the dynamics of a chain as a diffusion 
process [8]. The GNLuM may be applied to the cases of an isolated cold 
chain or of a hot polymer in the vapor phase. Applications of the GNLcrM 
have been developed in Refs. [HI [TO] . 

This work is organized as follows. In Sec.[2]the dynamics of a classical chain 
is investigated in three dimensions. The kinetic energy of a discrete chain 
with N — l segments is derived in cartesian and spherical coordinates. More- 
over, the limit to a continuous chain is performed. In Sec. [3] the probability 
distribution function for an inextensible chain in a heat bath is constructed 
using a path integral approach. Sec. [4] is dedicated to the discussion of the 
dynamics of a rigid chain in which the segments are allowed to form only 
fixed angles with respect to the z axis. Finally our Conclusions are drawn 
in Sec. (pjl). 

2. Classical dynamics of a three-dimensional chain with rigid 

constraints 

Let us consider a chain of N — 1 segments PjPj_i of fixed lengths li 
(i = 2, . . . , N) embedded in a three-dimensional space. With the symbol l\ 
we denote the distance of the end point Pi from the origin of the coordinate 
system. Additionally, there are small beads of mass mi attached at the 
joints of the segments PjPj_i, where i= 1, . . . , N. 

The above construction describes a freely jointed random chain, which is 
one of the basic models used in polymer physics. Freely jointed means that a 
given segment can take with equal probability any spatial orientation inde- 
pendently of the orientations of the neighbouring segments. The position of 
each segment PjPj_i can be specified by giving the coordinates of its endings 
Pi and Pi— i in cartesian coordinates Pi(t) = [%i{t),yi(t),Zi(t)]. However, in 
the following it will be more convenient to use spherical coordinates: 

i 

Xi(t) = lj cos (pj (t) sin 6j (t) (i = l,...,N) 
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i 

yi(t) = ^2 l j sirup j (t) sin 9 j(t) (i = l,...,N) 

3=1 

i 

= E Z J cos %(*) (i = l,...,JV) (1) 

3=1 

We will also neglect analytical complications connected with the inclusion 
of interactions such as the hydro dynamic interaction and steric effects. In 
this sense the chain is treated as a free one. 

The dynamics of a such a chain can be regarded as the motion of a 
system of coupled pendulums. For the sake of simplicity one of the ends 
of the chain has been fixed in the origin, see Fig. [T] Apart from that, no 
restrictions will be imposed on its motion. This implies that different parts 
of the chain are allowed to penetrate one into the other. In this case the 
chain is called a phantom chain. 

The fact that the chain is attached at the origin of the coordinates 
corresponds to the condition P\ = (0, 0, 0) or, equivalently: l\ = l\ = 0. 




Fig. 1. Three-dimensional conformation of a chain made of links and beads. 
The calculation of the kinetic energy of the system K^iscif) m spherical 
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coordinates is long but straightforward and gives as an upshot: 

N n-l m N n-l 

Kdisc = X! X ~2~ In-k+lfin-k+l sln ^n-k+l + X/ X/ ~^^n-k+l ®n-k+ 1 
n=l fc=l n=l fc=l 

A 1 " n— 1 n— A; 

+ X] X] X m nkln-k+l [<Pi<Pn-k+l sin 0* sin 6> n _ fc+ i cos (<£n-fc+l - <Pi) 
n=l k=l i=2 

+ 9itp n -k+i cos 6*i sin n _fc + i sin - v^^-fc+i) 
+ ^j^n-fe+i sin 0, cos n -k+l sin (<p n -fc+l - </?i) 

+ OiOn-k+i (cos 9i cos 6 n ^ k+ i cos (ip n -k+i ~<Pi) + sin 9i sin n ^ k+1 ) (2) 

To pass to the limit of a continuous chain we will use the rigorous procedure 
described in [7 a , where the two-dimensional case was analyzed. In order to 
do that, we assume that all segments have the same length and all beads 
have the same masses: 

M 

k = a (i = 2,...,N) m i = —a (i = l,...,N) (3) 

Li 

N 

where M = m i an d L are the total mass and the total length of the 

chain respectively. The next step consists in performing the limit in which 
the continuous system is recovered: 

a — ► N — ► +oo Na = L (4) 

One can see from Q that the product Na is fixed and gives the total length 
of the chain. Exploiting Eqs. (|2jj4| it is possible to get the kinetic energy of 
the continuous chain. Let's see for example how the recipe for performing 
the continuous limit works in the case of the third term in ([2]): 

N n—1 n—k 

XXX mnkln-k+lftiVn-k+i sin 6% sin n _ fe+ i cos(^ n _ fc+1 - tpi) 

n=l k=l i=2 

N n—1 n—k 

— > "F X a X a X ac ^' s n~ Sk+ a) 

n=l k=l i=2 

x sin 9(t, Si) sin8(t, Si — Sk + a) cos(cp(t, Si — Sk + a) — tp(t, si)) 
a —,Q M f L f s 

a = — ds(L — s) / dvip(t,v)(p(t, s) sin 9 (t,v) sin 6 (t, s) 
L Jo Jo 

x cos((p(t, s) — (p(t, v)) (5) 



directedchain printed on May 22, 2009 



5 



In obtaining the last part of the above equation we have exploited the for- 
mula 

/ ds I* duf(u) = [ ds{L-s)f{s) (6) 
Jo Jo Jo 

which is valid for any integrable function /(s). Applying the prescription of 
Eq. ([5j to the rest of the terms in Eq. ([2j , we get with the additional help 
of Eq. ([6| the full expression of the kinetic energy of the continuos chain: 



M 



L 



L ps 



JO 



K(t) = — ds{L — s) I dv (p(t,v)cp(t,s)sm6(t,v)sm0(t,s 



x cos ((f(t, s) — lf(t, v)) 

+ 9(t, v)(p(t, s) cos 9(t, v) sin 9(t, s) sin ((f(t, v) — (p(t, s)) 

+ 9(t, s)(p(t, v) sin 0(t, v) cos 9(t, s) sin (ip(t, s) — (p(t, v)) 

+ 0(t, v)9(t, s) ( cos 9(t, v) cos 9(t, s) cos ((p(t, s) — (p(t, v)) 

+ sin 9 (t,v) sin 0(t,s))} (7) 



We would like to stress that the right hand side of the derived equation 
contains five terms, while the initial discrete formula of the kinetic energy of 
Eq. ([2]) contained seven terms. This is due to the fact that the contributions 
from Eq. (jij) in which and 9^_ k+1 are present disappear after taking 

the continuous limit because they are proportional to a — ► 0. 

For future convenience, we give also the expression of the kinetic energy 
in cartesian coordinates: 

TV 

K disc = ^(x* + yl + z!) (8) 

i=2 Z 

where Xi, y% and Zi have been defined in Eq. ([T|. The sum over i starts from 
2 because one end of the chain coincides with the origin of the axes, so that 
l\ = 0. Of course, due to the condition that each segment has a fixed length 
li, Eq. Q must be completed by the following constraints: 

( Xi - Xi-x) 2 + ( yi - yi-xf + {zi - m-!) 2 = l\ (i = 2,...,N) (9) 

At this point we have thus two choices. Either we keep the kinetic 
energy in the simple form of Eq. ([8| at the price of having to deal with the 
constraints ^ , or we solve those constraints using the spherical coordinates 
li,9{, ifi of Eq. ([!]). In the latter case, the kinetic energy of Eq. is both 
nonlocal and nonlinear and thus is difficult to be treated. In the continuous 
limit, the situation does not change substantially. 
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We end up this Section performing the continuous limit of the kinetic 
energy of Eq. ([8]) and of the constraints ([9| . Following the prescriptions 
given in Eqs. pfPl), we obtain: 



M 
2L 



and 



o 



R /2 = l 

where we have introduced the vector notation: 

R = [x(t,s),y(t,s),z(t,s)] 



dsil 2 (t,s) 



(10) 

(11) 
(12) 



to describe the position on the chain. In polymer physics R is called the 
bond vector. In Eq. (10) and Eq. (El]) we have put R = ^ and R' = 
The compatibility of the description in cartesian coordinates with that in 
spherical coordinates can be verified by introducing the fields 9(t, s), </?(t, s) 
connected with the cartesian fields x(t, s),y(t, s), z(t, s) by the relations 



x(t, s) 
z(t,s) 



du cos ip(t, u) sin 9(t, u) 
du sin (p(t,u) sin 6(t, u) 
du cos 8(t, u) 



(13) 
(14) 
(15) 



If one performs the substitutions of Eqs. (13-15) in the kinetic energy (10) 
and makes use of the formula ( [6]) , one arrives exactly at the expression of the 
kinetic energy ([7]). Thus, Eq. (j7|) and Eq. (10) together with the constraint 
(11) are equivalent. 



3. Dynamics of a chain immersed in a heat bath 

In this Section the path integral formulation of an inextensible chain in 
the contact with a heat reservoir at temperature T is provided. According 
to the construction presented in Sec. [2] the conformation of the chain is 
treated as the limit case of a system of N beads connected by N — 1 links 
of fixed length a. In the discrete case the positions of the beads are given 
by a set of three-dimensional cartesian vectors R n (i), (rt = 2, . . . , N), while 
the conformation of the continuous chain at a given instant t is described 
by the vector field R(i, s), s being the arc-length. Furthermore, the chain 
is inextensible and thus has constant length L = Na. 

In order to describe the thermodynamic fluctuations of the chain, we 
regard it as a system of iV Brownian particles of mass m whose trajectories 
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satisfy the constraints of Eq. ([9]). These constraints enforce the condition 
that the total length of the links connecting the beads should be equal to 
a. It is possible to rewrite Eq. ([9| in the more compact form: 

|R„(i) - R n _i(i)| 2 

5 = 1 (n = 2,...,N) (16) 

cr 

We also require that at the initial and final times t = and t = tf the 
position of n-th particle is respectively given by R n (0) = Ro,n and TL n (tf) = 
R/, n for n = 2, . . . , N. 

In other words, the primary task of this Section is to analyze the dynam- 
ics of a system which consists in the constrained random walk of the beads 
composing the chain. The main difficulty in performing analytical calcula- 
tions are obviously the constraints. Starting like in the Rouse model from 
an approach to the problem based on the Langevin equation to describe 
the motion of a polymer in a solution |llj , the treatment of the constraints 
becomes awkward. For this reason we will use an interdisciplinary strat- 
egy, which combines the techniques of field theory with those used in the 
statistical mechanics of polymers with topological constraints. The starting 
point of the presented framework is to specify the probability distribution 
function expressed in a path integral form. contains the physical 
information about the system. To be more specific, it measures the proba- 
bility that the chain after a given time t / passes from an initial configuration 
Ro, n to a final configuration R/, n . 

Before we construct the probability function for the chain with rigid 
constraints, let's see how the path integral of a single free Brownian particle 
looks like. In order to do this we assume that at the time t = the particle 
finds itself at the initial point Ro and starts to perform a random walk. As 
it is well known, the probability V(^/; R/> Ro) that, after the time tf the 
particle arrives at a given point Rj, satisfies the diffusion equation 

dt f 9R 2 { ' 

where D is the diffusion constant. The boundary condition at tf = is 
chosen in such a way that ^(0,R/;Ro) = <5(R/ — Ro)- The solution ip of 



(17) can be expressed in the form of a path integral 

1 R ' 1W , /-*/R 2 (t) 
PR(t)exp - / 

R(0)=R 



V>(*/,R/;Ro) = A / DR(t)exp -/ 

JR(0)=R JO 



dt 



(18) 



AD 

where A is a normalization factor. We note that the diffusion constant D 



appearing in Eq. (18) satisfies the relation D = where k is the Boltz- 
mann constant, T is the temperature of heat bath and r is the relaxation 
time that characterizes the rate of decay of the drift velocity of the particle. 
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The above prescription can easily be generalized to a system of N non- 
interacting Brownian particles. It this case the probability that the n— th 
particle starting from the point Ro, n arrives at the point R/, n is given by 



N 



N 

n 

71=1 



R„(t/)-R;,„ 

A I T>Il n (t) 

R„(0)=R O ,„ 



exp 



1 



2kTr 



N 

E 

71=1 



tf 



(19) 



In addition, the form of the path integral on the right hand side of Eq. ( 19 ) 



displays the connection with the partition function of a set of N free par 

N 

tides in quantum mechanics where the functional An — z-, jo ~2 

n=l 

represents the action of the system. The well known duality between quan- 
tum mechanics and Brownian motions allows to treat the factor 



E So' ^n(t)dt 



K = 2/cTr 



(20) 



as the quantity which plays the role of the Planck's constant. Indeed, one 
may show that the uncertainties in the position and momentum of a Brown- 
ian particle due to the frequent collisions with the molecules in the solutions 
satisfy an analog of the Heisenberg uncertainty relations: ApAr ~ k [12 . 

Going back to the dynamics of an inextensible chain, the only difference 
with respect to a system of free particles is that the bond vector R n (t) 
satisfies the additional constraints (16) restricting the trajectories of motion. 



To implement them in the dynamics of noninteracting Brownian particles 



we add a product of functional delta functions in the path integral (19) 
which imposes the desired conditions (16): 



N 



c 

N 



n / pR -w 

_^=l JRn(0) = R ,„ 

2 



X W |Rn(*)-Rn-l(*)|' _ 1 

n=2 V a / 



(21) 



where C is an irrelevant factor and the mass of a single particle present in 



Eq. (19) has been replaced according to the equation m = ^f-a. The above 



procedure to fix the constraints in a path integral has been applied in the 

statistical mechanics of entangled polymers |13| [T^l [15] . 

The next step consists in performing the continuous limit Q in (21) |16j : 



R(t f ,s)=R f (s) 
R(0,s)=R (s) 



PR(t,s)e"5sWi'o /dt s/o id ^ 2 (^)5(R /2 (t, s ) - 1) (22) 



The result of Eq. (22) defines a model which is closely related to the non- 



linear sigma model (NLerM) used in high energy physics [T7], solid state 
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physics [18] and disordered systems [19]. For this reason it has been called 
generalized nonlinear sigma model (GNLcrM). The most striking difference 
between these two models lays in the constraints, which in the case of the 
NLcjM are of the form R 2 = 1, while in the GNLcrM they have been replaced 
by the nonholonomic condition (11). 

To conclude this section let us note that it is possible to show that the 
generating functional of the correlation functions of the GNLcrM coincides 
with the generating functional of the correlation functions of the solutions 
of a constrained Langevin equation |8 . 



4. Dynamics of an inextensible chain with constant bending angle 

The approach presented in Sec. [2] in order to treat the dynamics of 
random chains has some interesting variants which we would like to discuss 
in this Section. To this purpose, we choose the formulation in which the 
positions of the ends of the segments composing the chain are given in 
cartesian coordinates. As we have already seen, in this way the expression 
of the kinetic energy K^ isc is given by Q and must be completed by the 
constraints ([9]). From now on we assume as before that all segments have 
the same fixed length l n = a, but we require additionally that: 

(z n - zn-x) 2 = b 2 <a 2 (23) 

This implies that the projection of each segment onto the z— axis has length 
±6, so that the segments are bound to form with the z— axis the fixed angles 
ai = a or «2 = (vr — a) defined by the relations: 

b b 
cosai = +- cosct2 = — (24) 

a a 



Clearly, in both cases the constraints (|9j) and (23) may be rewritten as 
follows: 

( Xn -Xn-l) 2 + (jjn ~ Vn-l? = _J_ _ 1 ( n = 2,...,iV) (25) 

b z b z cos z a 

where a may be either a\ or 012 ■ In the following we will suppose that 
only the angle a = ol\ is allowed, so that the chain cannot make turns in 
the z direction. An example of a conformation of a chain satisfying these 
assumptions is given in Fig. [2] 

The constraints (25) are eliminated using the spherical coordinates of 
Eq. ([!]) after setting the angles 9 n formed by the segments with the z— axis 
equal to a: 



n 



x n{t) = 2_. cos Lpi{i) sin ol (26) 



i=l 
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Fig. 2. Example of motion of a chain whose segments are constrained to form a 
fixed angle a with the z— axis. In the figure a = 30°. 



Unit) = ^2 k sin ipi(t) sin a 



i=l 
n 



= cos a = na cos a 



(27) 
(28) 



i=i 



As we see from the above equation, each segment is left only with the 
freedom of rotations around the z— direction, corresponding to the angles 
tpi(t). Moreover, the total length of the chain is always L = Na, but now 
also the total height h of the trajectory along the z— axis is fixed: 



h = Nb 



(29) 



At this point, we pass to the continuous limit, this time taking as parameter 
describing the trajectory of the chain the variable z instead of the arc-length 



s. Due to the last of Eqs. (28), the z— components of the velocities are always 
zero: 

z n (t) = (30) 

As a consequence, we are left with something similar to a two-dimensional 
problem. The difference from a real two-dimensional problem, which could 
be obtained by putting 9j = -ir/2 (j = 1,...,N) in Eq. ([!]), is that the 
equations describing the position of a bead in two dimensions, namely 
Xi(t) = Yjf=i h CQS <£i(t) an d Ui(t) = YljLi h shxfj{t), have been replaced 
by Eqs. (26) and (27). Moreover, the constraints have a slightly different 
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form. Following the same procedure presented in Sec.[2j we find after a few 
calculations the expression of the kinetic energy in the continuous limit: 

M f h f z f Zl 
K a = tan a— / dz dz\ \ dz-i 
2h Jo Jo Jo 

x tp(t, z — zi)<p(t, Z2) cos(tp(t, z — zi) — (f(t, Z2)) (31) 



and of the constraint (25 ): 

(3 z x) 2 + (d z y) 2 = tan 2 a (32) 

It is also not difficult to show that the probability distribution \& Q is given 
in cartesian coordinates by: 



= J Vx{t, z)Vy{t, z) exp{-^}<5((d z x) 2 + (d*y) 2 -tan 5 



a) (33) 



x 2 (t,z)+y A (t,z) (34) 



where 

M f l f , f h , 
A* = ttt / dt dz 
2ai Jo Jo 

At this point we discuss briefly the case in which both angles tt — a and 
a are allowed. In this situation, the trajectory of the chain may have turns. 
An example of motion of this kind is given in Fig. [3} 

The constraints ^ and (23) remain unchanged, but the coordinate z 



cannot be chosen as a valid parameter of the trajectory of the chain and one 
has to come back to the arc-length s. The most serious problem is the fact 
that the variables z n {t) are not continuous functions of the time, since the 
length z n (t) — z n -\{t) is allowed to jump discretely between the two discrete 
values +b and —b, corresponding to the angles a and tt — a respectively. It 
is therefore difficult to define the components z n of the velocities of the ends 
of the segments and thus their contribution to the kinetic energy. Let us 
note that this problem affects only the z degrees of freedom. The degrees 
of freedom x n {t) and y n (t) of the chain remain continuous functions of t 
despite the jumps of the z n 's. This fact can be easily verified looking at the 



definition of x n (t) and y n (t) in Eqs. (26) and (27). Since sin(7r — a) = sina 



both the x n (tys and y n (i)'s are not affected by the jumps of the angle 
a < — * tt — a. As a consequence, the problems with the z variable can easily 
be solved if the chain has no interactions in which the z variable is involved. 
In this case, in fact, the degrees of freedom connected to the motion along 
the z— direction are decoupled from the other degrees of freedom and may 
be neglected. 

As a consequence, we assume that the interactions are z— independent, 
so that the difficulties related to the motion along the z— direction disappear 
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Fig. 3. Example of motion of a chain whose segments are constrained to form 
fixed angles a or n — a with the z— axis. In the figure a = 30°. Turning points are 
emphasized by means of shaded beads. 



and once again the problem reduces to that that of a two-dimensional chain. 
Since the constraints are always those of Eqs. ^ and (23 ), one may proceed 
as in the case of fixed angle a. As a result, one finds that the final probability 
distribution is of the form: 



C J Vx(t, s)T>y(t, s) exp 



where 



Aa 7T- 



M 
2L 



dt 



5{{d s x) 2 + {d s y) 2 - tan 2 a) 



ds 



x 2 (t,s)+y 2 (t,s) 



(35) 



(36) 



and C is a constant containing the result of the integration over the decou- 
pled z degrees of freedom. With respect to the previous case, let us note 



that in Eqs. (35) and (36) z has been replaced by the arc-length s as the 



parameter of the trajectory of the chain. Correspondingly, the total chain 
length L appears in the action instead of the height h. 



5. Conclusions 



In this work the dynamics of an inextensible freely jointed chain con- 
sisting of links and beads in three dimensions has been discussed both from 
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the classical and statistical point of view. In Sec. [2] we have mainly concen- 
trated ourselves on the computation of the kinetic energy in cartesian and 
spherical coordinates for the discrete and continuous chain. In Sec. [3] we 
have derived the probability function of the chain using a path integral 
framework and the fact that the fluctuations of the chain can be regarded 
as those of a system of Brownian particles with an additional constraint 
condition imposed on their trajectories. The probability function ^ of this 
system is equivalent to the partition function of a generalized nonlinear a 
model. The analogy of the GNLoM with the NLcrM suggests the possibility 
of applying techniques and results coming from the NLcM to the GNLuM. 
For example, it is known that the NLcrM is renormalizable in two dimen- 
sions and also that it has interesting features because it is analytically free 
and has a dynamically generated mass gap [20]. The similarity with NLcrM 
seems also to suggest that there is no symmetry breaking in the underlying 
0{d) symmetry of the GNLcrM, where d denotes the dimension of the vector 
field R(i, s). One should however be careful when extending the results of 
the NLcrM to the GNLuM. For example if d = 2, one may use polar field 
coordinates to express the vector field R(i,s). If one does that the NLcrM 
becomes a free field theory in the angle variable |21J . This is not true in 
the case of the GNLuM which in polar coordinates exhibits a nonlinear and 
complicated dependence on the angle degree of freedom. Moreover, it is 
not straightforward to apply techniques like the effective potential method 
which is useful to investigate possible phase transitions in the NLcrM. The 
reason is that in this method it is performed an expansion around field con- 
figurations minimalizing the action which are constant. Configurations of 
this kind correspond in the GNLuM to the situation in which the chain has 
collapsed to a point and thus are nonphysical. Finally in Sec. [4] a three- 
dimensional chain admitting only fixed angles with respect to the z— axis 
has been discussed. It has been shown that it is possible to reduce the 
problem to two dimensions, in a way which is similar to the reduction of 
the statistical mechanics of a directed polymer to the random walk of a 
two-dimensional particle [22 . Our approach is valid only if the chain has 
no turning points. If there are turning points the kinetic energy is not well 
defined, because the variable z(t, s) is no longer a continuous function and 
thus its time derivative becomes a distribution. One way for adding to our 
treatment turning points as those of Fig. [3] is to replace the variable z with 
a stochastic variable which is allowed to take only discrete values. Another 
way is to look at turning points as points in which the chain bounces against 
an invisible obstacle. A field theory describing a one-dimensional chain with 
such kind of constraints has been already derived in Refs. [23] . 
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